Propagating wave-packets and quantised currents in coherently driven polariton 

superfluids 
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We study the properties of propagating polariton wave-packets and their connection to the sta- 
bility of doubly charged vortices. Wave-packet propagation and related photoluminescence spec- 
tra exhibit a rich behaviour dependent on the excitation regime. We show that, because of the 
non-quadratic polariton dispersion, doubly charged vortices are stable only when initiated in wave- 
packets propagating at small velocities. Vortices propagating at larger velocities, or those imprinted 
directly into the polariton optical parametric oscillator (OPO) signal and idler are always unstable 
to splitting. 

PACS numbers: 42.65.Yj, 47.32. C-, 71.36.-l-c 



After almost a decade since the first reports of bosonic 
stimulation [l|, y], resonantly pumped polaritons have 
been shown to exhibit a new form of non-equilibrium 
supcrfluid behaviour [3|-|5|- The experimental ability to 
control externally the polariton currents (phase varia- 
tions) with the laser pump [J] as well as to generate 
propagating polariton wave-packets (density and phase 
variations) , by applying another pulsed laser probe [3| , is 
a unique feature of driven polariton systems. Here, the 
interplay between phase and amplitude variations opens 
the possibility for novel superfluid phenomena. In ad- 
dition, by using a pulsed Laguerre-Gauss beam, polari- 
tons in the optical parametric oscillator (OPO) regime 
have been recently shown to display persistence of cur- 
rents and metastability of quantum vortices with charge 
171=1 and 771 = 2 [5|. 

In this Letter we report a comprehensive theoretical 
analysis of propagating polariton wave-packets, triggered 
by a pulsed probe, and examine the stability of vortices 
of charge m ~ 2 in different excitation regimes. In par- 
ticular, we asses how the wave-packet shape, velocity and 
intensity, as well as the stability of to = 2 vortices, are 
related to the probe wavevector and to the photolumi- 
nescence (PL) spectra. The behaviour of to = 2 vortices 
has been studied recently Q, and two distinct regimes 
have been singled out: In the TOPO regime the initiated 
vortex only lasts for as long as the additional decaying 
gain introduced by the triggering probe. Here, to = 2 
vortices have been shown to be remarkably stable when 
initiated at small group velocities and instead to split 
into two TO = 1 vortices at larger velocities. In a second 
regime, vortices have been shown to withstand the gain 
and imprint into the OPO. In this case, m = 2 vortices 
always split. We provide the theoretical explanation for 
this intriguing phenomena. 

The behaviour of ttt, = 2 vortices has been the subject 
of intensive research in the context of ultra-cold atomic 



gases. However, although stable free to = 2 vortices 
have been predicted for specific ranges of density and 
interaction streiigth [6|, |7[ , they have not been observed 
experimentally [8|. Stable pinned to = 2 persistent cur- 
rents have been recently realised [9| only by using toroidal 
traps. 

Model The time evolution of coherently driven po- 
laritons is described by the Gross-Pitaevskii equations 
for coupled cavity photon and exciton fields il)c,x{'<^,t) 
with pump and decay {h = 1) [lO|: 
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We use here the same notation as in Rcf. jll|. The cav- 
ity field is driven by a continuous- wave pump. -Fp(r, t) = 
•^i ,c {'i')e^'^^f ' ""-'^p*) with a Gaussian or a top-hat spa- 
tial profile J-f ,cr of strength fp and full width at half 
maximum (FWHM) Up. Fpb{v,t) is the pulsed probe in- 
troduced later. We neglect the exciton dispersion and as- 
sume a quadratic dispersion for photons, ujc ~ ^[7^2 — ■ 
rifl is the Rabi splitting and the fields decay with rates 
Kx,c- The exciton interaction strength gx is set to one by 
rescaling both fields ipx,c and pump strength Fp and Fph 
by y/Qii/{2gx)- Through the paper, mc = 2 x 10~^too, 
the energy zero is set to ujx = w^ (zero detuning) and 
il.R ~ 4.4meV. We numerically solve Eq. ([T|) on a 2D grid 
using 5"^-order adaptive-step Runge-Kutta algorithm. 

In the absence of the probe {Fpt = 0), polaritons are 
continuously injected into the pump state and, above the 
pump strength threshold for OPO, fp^, undergo stim- 
ulated scattering into the signal and idler states [^ 
(see Fig. [T] left panels). We filter the signal, idler and 

pump spatial profiles |^^*^(r, i)|e"*'<^-^ ^'"'*' in a cone 
around their associated momenta i^s,i,p- The supercur- 
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FIG. 1: (Color online) PL spectra (top) and spatial profiles of pump and filtered signal 1^^0(^10,^)1 (bottom) for three regimes: 
(i) OPO (left panels) for a smoothen top-hat pump above threshold (fp = 1.12/p ): polaritons, continuously injected at 
{kp,0) and ijp, undergo stimulated scattering into signal and idler stationary steady states {kx,c = 0.2meV). (ii) TOPO 
(middle panels): a short, at = Ips, Laguerre-Gauss m = 2 (top) and Gaussian m = Q (bottom) probe ([5} resonant with (fcp(,,Q) 
and u)pb triggers propagating signal and conjugate states, which lock to the same group velocity (kx = 0, kc = 0.02meV). 
(iii) Similar conditions to (ii) but k.x,c = 0.09meV (right panels); now there is no amplification of the probe, both signal and 
conjugate decay exponentially. Lower (LP) and upper (UP) polariton (blue dashed lines), and bare cavity photon (C) and 
exciton (X) dispersions (gray dashed lines) are shown for reference. 



rent, V(/)^*^, is a superposition of the dominant uniform 
flow \<is.i,p a-nd more complex currents (caused by the sys- 
tem being finite size), which move particles from gain 
to loss dominated regions. In addition one can charac- 
terise the density variations defining the group velocity 
^s,'i,p _ dj-^.hp ^di^ where rf^'P is the maximum of the 
signal (idler and pump) spatial profile. The OPO is a 
steady state regime, where \ip^^^{r,t)\ are time indepen- 
dent. This is also recognisible in the typical flat disper- 
sion of signal, idler and pump PL spectra (see left panel 
of Fig. [1]) — the signal (idler, pump) group velocity being 
the derivative of the dispersion at k,, (kj;,p). 

Propagating wave-packets and TOPO regime In order 
to create propagating wave-packets (finite group veloc- 
ity) , one needs to add a pulsed probe with finite momen- 
tum kpb, spatially small with respect to the pump spot 
size [3|. Since we are interested in the stability of trig- 
gered vortices of charge m = 2, and its connection to the 
character of the propagating wave-packet, we consider 
here a Laguerre-Gauss pulsed probe (sl. Illl. Il2| . 
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where the phase ip winds from to 27r around the vortex 
core Tpb- Above threshold for OPO, fp > fp^, one can 



identify two distinct regimes depending on the system pa- 
rameters |5|: In one case, the triggered vortex propagates 
as extra population (gain) on top of the OPO stationary 
signal and idler states [3[ — TOPO regime. Here, the 
vortex lasts only for as long as the additional decaying 
gain. A similar TOPO regime can be found also below 
the OPO threshold, and for simplicity we show this case 
in Fig. [1] (middle panel). Alternatively, vortices persist 
the gain and get imprinted into the OPO signal and idler 
— metastable vortex regime [11| . 

In the TOPO regime, the polaritons generated by the 
probe scatter parametrically with the ones of the pump, 
generating a traveling signal (here meant as extra popula- 
tion on top of either the OPO signal or idler, or full signal 
if below the OPO threshold) with momentum k^ = kpf, 
(which can be either kpj, > kp or kpf, < kp) and its conju- 
gate at kc — 2kp — kpf,. In the case of the middle panels of 
Fig. [T] the signal and the conjugate get initially strongly 
amplified by the parametric scattering from the pump 
state, then decay slowly (we refer to this as a proper 
TOPO regime) and, at later times, eventually the decay 
becomes exponential. This regime reproduces the one ob- 
served in experiments (see Fig. 3 of Ref. [14|), where the 
intensity maximum is reached quickly, within 4ps after 
the maximum of the pulsed probe, and is immediately 
followed by a slow decay. Instead, for the case of the 
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FIG. 2: (Color online) Left panel: Group velocity, v"g of the 
propagating signal wave-packet as a function of the probe 
momentum fcpj,. The black dots are determined from simula- 
tions, whereas the solid line (red) is the derivative of the LP 
dispersion evaluated at fcpf,, v]^\- The green dashed line in- 
dicates where the LP dispersion deviates from the quadratic 
behaviour. Right panel: Signal (solid lines) and conjugate 
(dashed lines) states for fcp(, = 1.4/im^^ at different times 
after the arrival of the probe. 



right panels of Fig. [TJ parametric scattering is too weak 
to induce any significant amplification, and we observe 
exponential decay of signal and conjugate populations 
immediately after the probe switches off. 

By analysing the change in time of the spatial pro- 
file of the TOPO signal |i/;J^(r, i)|, wc find that its 
group velocity v* is simply given by the derivative of 
the lower polariton (LP) dispersion evaluated at kph, 
i.e. for zero detuning and low densities by ii^^ = 

/cpfc/(2mc) - /sV(2mcV^'' +4m^0^) (see Fig. [5]) . We 
have confirmed this for both TOPO regimes. 

In order to understand this result it is instructive to 
analyse the PL spectra. In the regime with weak para- 
metric scattering (right panel of Fig. [T|), aside the strong 
emission from the pump state, the dispersion is simply 
that of the LP, and thus the signal propagates with a 
group velocity given by v^^ . However, in the proper 
TOPO regime (middle panel of Fig. [T]) the spectrum be- 
comes linear. This can be explained by the fact that in 
order to have efficient parametric scattering processes, 
signal and conjugate require a large spatial overlap and 
so their group velocities need to lock to the same value, 
therefore the dispersion becomes linear — a similar result 
has been found in ID simulations, as well as in experi- 
ments [3[. Here, we determine that the TOPO linear 
dispersion is tangential to the LP branch at fcp;,, thus its 



slope is given in this case also by w^^. 

The linearisation of the TOPO dispersion, with the 
slope defined by the probe momentum fcpb, also explains 
our observation that the TOPO amplification is stronger 
when signal (k^ = kpf,) and conjugate (kc) wavevectors 
(connected via parametric constraint 2kp = kc,-|-kc) have 
nearby values. This implies that the LP tangential line 
at kpb lies close to the LP dispersion also at kc, therefore 
leading to efiicicnt stimulation. It is interesting to note, 
that by changing kp and kp(, it is possible to engineer 



wave-packets with a broad range of values of group ve- 
locities and supercurrents — e.g., for kpb = 2kp, one has 
a conjugate state at kc = (i.e., no net current) with a 
large group velocity, u^ = v^^ . 

From the PL spectrum we can deduce the nature of 
the wave-packet propagation. For linear dispersion one 
expects a soliton-like behaviour, with signal and conju- 
gate propagating and not changing neither shape nor in- 
tensity. For quadratic dispersion, propagation is anal- 
ogous to the ballistic time-of-fiight expansion of ultra- 
cold atomic gases [l3J, i.e. for a Gaussian wave-packet, 
the FWHM=(ct2^ -F ( zmc^ y 'f^ grows in time, but the 
shape is preserved. The total density decays exponen- 
tially with rate given by '^'^^^ at zero detuning, and 
the maximum of the wave-packet moves with a constant 
velocity v\^^^ (equal to -^^^ in the quadratic LP regime). 
Due to the dynamical nature of the TOPO state, the sys- 
tem evolves between different regimes. In particular, only 
in the strong amplification regime the spectrum is linear, 
while it evolves back to LP at longer times. In addition, 
when the LP dispersion deviates from quadratic, propa- 
gation becomes complex: the wave-packet gets distorted 
and one observes beatings in the spatial profiles. 

Propagation of signal and conjugate are shown in the 
right panel of Fig. [51 We observe a mechanism analogous 
to the one found in four-wave-mixing experiments |15l . 
|l6|: when the probe arrives the conjugate propagates 
faster then the signal, before getting locked to it with a 
small spatial shift of their maximum intensities. At later 
times, when the density drops and the parametric process 
becomes inefficient, the two wave-packets start unlocking 
— the conjugate slows down with respect to the signal if 
fcc < kg as in Fig. [21 or moves faster when kc > ks- 

Stability o/ m = 2 TOPO vortices In the TOPO 
regime, we observe m = 2 vortices to be stable within 
their lifetime when triggered at small momenta kpb 
(Fig. [3K,c), while TOPO vortices always split into two 
TO = 1 vortices for large kpb (Fig. [Sbjd), in agree- 
ment with the recent experimental observations [5|. Our 
numerical analysis shows that the crossover from non- 
splitting to splitting happens for values of the probe mo- 
mentum where the LP dispersion deviates from quadratic 
(see Fig. [2]). The two regimes are shown in Fig. [3l For 
^ph=0.2^m~^, at short times, the signal carrying the vor- 
tex propagates without changing shape and with little 
change in intensity, consistent with the corresponding 
PL spectrum having a linear dispersion. However, at 
longer times the density drops more then two orders of 
magnitude, the dispersion becomes quadratic and there- 
fore the wave-packet expands (panel c). A uniform ex- 
pansion of the wave-packet leads to the increase of the 
vortex core but does not cause the vortex to split. In 
contrast, for fcpf,=1.4^m~^, where the LP dispersion is 
not quadratic, the splitting of the m = 2 vortex state 
into two TO = 1 vortices happens shortly after the arrival 
of the probe (panel d) . This behaviour can be explained 



100 




-100 



t = 150ps 



100 

X [)im] 



r # V 



-100 

X [[im] 




-10 10 -10 10 

X [^m] X [jjim] 

FIG. 4: (Color online) Filtered signal profile and currents 
above OPO threshold fp = 1.12/p for a top-hat pump with 
FWHM ap — 35/im at i = —3ps (left panel) before the arrival 
of an m, = 2 probe at f = Ops (right). The double quantised 
vortex splits into two m = 1 vortices even before the probe 
reaches it's maximum intensity. 



FIG. 3: Profiles of TOPO signal intensity |7/>J(r, t)\ and phase 
</>cr(r, t) (inset) after the arrival (at t = 0) of an m = 2 vortex 
pulsed probe ([2} with apb — 87^m. The conditions are the 
same as the middle panels of Fig. [T] but with kpt = 0.2^m~^ 
(a,c) and kpb = 1.4/xm~^ (b,d). While in the case (a,c) the 
m — 2 vortex does not split within its lifetime, in (b,d) the 
vortex splits soon after the probe arrives. The intensity scale 
in (c) is 200 times smaller then in (a) - signal expands as the 
density drops two orders of magnitude (see text). 



as follows: the dispersion of the time-dependent TOPO 
evolves from LP (before the probe arrival) to linear (at 
early times after the probe arrival), and back to the LP 
dispersion at later times. Wave-packets propagating with 
non-quadratic dispersion do not keep their shapes. In- 
deed, we see the distortion to be very pronounced at 
later times of the evolution. This distortion during the 
propagation leads to the mechanical splitting of m = 2 
vortex, analogous to the structural instability discussed 
in [17{ . Also, as discussed in [5|, for small kpb, within the 
quadratic part of the dispersion, the group velocity of 
the propagating vortex and the velocity associated with 
the net supercurrent (given by kpb) are equal. However, 
this is not the case for larger kpb, beyond the quadratic 
part, and so the propagating vortex feels rather strong 
net current in its frame, which may provide additional 
explanation for splitting. 

m = 2 vortices imprinted into the OPO As in Ref. [5| , 
we find that doubly charged (m = 2) vortices which get 
imprinted into the OPO signal are never stable, and split 
into two m = 1 vortices almost immediately, even before 
the probe reaches its maximum intensity, as shown in 
Fig. m There are several causes of the splitting: Before 
the probe arrives, the OPO dispersion is flat around the 
signal, idler, and pump. However, the triggering probe 



favours the signal and conjugate to lock and propagate 
with the same velocity v^^ , which corresponds to a lin- 
ear dispersion. Further, once the vortex gets imprinted 
into the stationary OPO signal and idler, the dispersion 
changes back to be flat. The evolution of the dispersion 
between linear and flat leads to a complicated dynamics 
of both signal and idler (the transient period described 
in [Ul), causing the structural instability and splitting 
of the vortex during the transient time. Another reason 
for the structural instability during the vortex imprinting 
into the OPO signal and idler are the non-uniform cur- 
rents (shown in Fig. 2]) caused by the interplay between 
spatial inhomogeneity, pump and decay, which the OPO 
vortex experiences in its reference frame. 

To conclude, we have determined group velocity and 
supercurrents of propagating polariton wave-packets trig- 
gered by a short pulsed laser probe, and established the 
conditions under which triggered m = 2 vortices are sta- 
ble and when they split into two vortices oi m = 1. Both 
phenomena are related to the form of the PL spectrum 
as well as the wave- vector of the triggering probe. The 
ability to control the topological charge and number of 
vortices, simply by changing the wave-vector of the ex- 
ternal probe, holds potential for operations in quantum 
information. In addition, we have shown that, by varying 
the pump and probe wavevectors, it is possible to engi- 
neer polariton wave-packets with any relation between 
group and supercurrents velocities — a property which 
can be useful for storage applications. 
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